

Dynamics of Sheared Granular Materials 

Lou Kondic 

Department of Mathematical Sciences and Center for Applied Mathematics & Statistics, 
New Jersey Institute of Technology, Newark, NJ 07102 

Brian Utter and Robert P. Behringer 

Department of Physics and Center for Nonlinear and Complex Systems 
Duke University, Durham, NC 27708 


ABSTRACT 

This work focuses on the properties of sheared granular materials near the jamming 
transition. The project currently involves two aspects. The first of these is an experiment 
that is a prototype for a planned ISS flight. The second is discrete element simulations 
(DES) that can give insight into the behavior one might expect in a reduced-g environment. 

The experimental arrangement consists of an annular channel that contains the granular 
material. One surface, say the upper surface, rotates so as to shear the material contained in 
the annulus. The lower surface controls the mean density/mean stress on the sample through 
an actuator or other control system. A novel feature under development is the ability to 
Thermalize 1 the layer-i.e. create a larger amount of random motion in the material, by using 
the actuating system to provide vibrations as well control the mean volume of the annulus. 
The stress states of the system are determined by transducers on the non-rotating wall. 
These measure both shear and normal components of the stress on different size scales. 
Here, the idea is to characterize the system as the density varies through values spanning 
dense almost solid to relatively mobile granular states. This transition regime encompasses 
the regime usually thought of as the glass transition, and/or the jamming transition. 

Motivation for this experiment springs from ideas of a granular glass transition, a related 
jamming transition, and from recent experiments. In particular, we note recent experiments 
carried out by our group to characterize this type of transition (Howell et al. Phys. Rev. Lett. 
82, 5241 (1999)) and also to demonstrate/characterize fluctuations in slowly sheared systems 
(Miller et al. Phys. Rev. Lett. 77 , 3110 (1996)). These experiments give key insights into 
what one might expect in near-zero g. In particular, they show that the compressibility 
of granular systems diverges at a transition or critical point. It is this divergence, coupled 
to gravity, that makes it extremely difficult if not impossible to characterize the transition 
region in an earth-bound experiment. 

In the DE modeling, we analyze dynamics of a sheared granular system in Couette 
geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better 
understand what we might encounter in a reduced-g environment, and at a deeper level to 
deduce the physics of sheared systems in a density regime that has not been addressed b\ 
past experiments or simulations. 

One aspect of the simulations addresses sheared 2D system in zero-g environment. For low 
volume fractions, the expected dynamics of this type of system is relatively well understood. 
However, as the volume fraction is increased, the system undergoes a phase transition, as 
explained above. The DES concentrate on the evolution of the system as the solid volume 
fraction is slowly increased, and in particular on the behavior of very dense systems. For 
these configurations, the simulations show that polydispersity of the sheared particles is a 
crucial factor that determines the system response. Figures 1 and 2 below, that present the 
total force on each grain, show that even relatively small (10 %) nonuniformity of the size 
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of the grains (expected in typical experiments) may lead to significant modifications of the 
system properties, such as velocity profiles, temperature, force propagation, and formation 
of shear bands. 

The simulations are extended in a few other directions, in order to provide additional in- 
sight to the experimental system analyzed above. In one direction, both gravity, and driving 
due to vibrations are included. These simulations allow for predictions on the driving regime 
that is required in the experiments in order to analyze the jamming transition. Further- 
more, direct comparison of experiments and DES will allow for verification of the modeling 
assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of 
these simulations of an annular system in zero-g environment will conclude the presentation. 



Figure 1: The (instantaneous) force experienced by the variable size particles; the volume 
fraction at the instant shown is approximately 0.9. The flow is driven by constant velocity 
motion of the upper wall in the +x direction (left to right). 



Figure 2: The force in the case of uniform size particles; all the other quantities are the same 
as in Fig. 1. Note formation of a shear band at y ~ —0.1; it should be mentioned that this 
particular shear band is just one of possible configurations characterizing the flow, and the 
manner in which the force propagates. 
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Main goal 

Understand the dynamics of a granular system 
as the volume fraction 7 is modified between 
gas-like and solid-like 

Analyze the phase transition which may occur 
in 3D zero-g experiment 

Analogy to 2D? - see Howell et al. PRL (1999) 
Discuss the parameters defining the system, in 
particular distribution of particle sizes 

Techniques 

Experimental 

Computational (discrete element simulations) 

This talk: Simulations 

Discuss expected behavior in zero g 
Analyze force propagation 
Also: can we thermalize the system in non- 
zero g and expect to reach uniform states? 
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S/(£2D/d) 




3900 particles 
volume fraction: 80% 



Experiment: 

2D: Second-order phase transition as the vol- 
ume fraction is decreased 

2D: gravitational effects are reduced by using 
a horizontal system 

3D (annular geometry): gravitational compaction 
does not allow to analyze dynamics close to 7,- 
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(an example of gravitational compaction) 

Goal of Simulations 

Analyze open (channel) (2D) Couette geom- 
etry to understand the basics of dynamics in 
3D experiment in annular geometry and zero -g 
environment 
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• Sheared system with slowly varying 7 

• Switch on gravity as well as excitation 
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From exponential to linear velocity profiles 
(note: linear profiles are not observed in annu- 
lar geometry!) 

From gas-like to solid-like behavior: jamming 
Significant differences between monodisperse 
and polydisperse systems: 

Fracture occurs for monodisperse systems 
Monodisperse systems are compressed easier 
(crystallization) — * large difference in the stored 
elastic energy 

Rate-dependent behavior for large volume frac- 
tions 
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temperature / vel. shear 2 elast. ener. / vel. shear 2 


Elastic energy 



Temperature 
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Sheared system + gravity 

Excitations through vibrated lower boundary 


Analyze: 
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Plain Couette Flow 

volume fraction: 80% (last 5 averages shown) 

variable size particles 
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Answer: It is not easy to fight Gravity! 

Difficult to produce uniform and dense granu- 
lar system which responds to shearing 
(there are also large forces resulting from ex- 
citations) 

Suggestion of this study: will need to perform 
experiments (and simulations) in zero g to fully 
understand the problem 


NAS A/CP— 2002-2 1 1 2 1 2/VOL I 


470 



Analyze in more detail the role of external exci- 
tations in both zero and nonzero g environment 
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Flow under gravity, shearing and excitations 


variable size particles 
Gamma = 4 
amplitude = d 
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